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Abstract 

We study radiative corrections to deeply virtual Compton scattering in the kinematics of HERA 
collider experiments to next-to-leading and next-to-next-to-leading order. In the latter case the 
radiative corrections are evaluated in a special scheme that allows us to employ the predictive 
power of conformal symmetry. As observed before, the size of next-to-leading order corrections 
strongly depends on the gluonic input, as gluons start to contribute at this order. Beyond next- 
to-leading order we find, in contrast, that the corrections for an input scale of few GeV 2 are 
small enough to justify the uses of perturbation theory. For £ > 5 ■ 10 -3 the modification of the 
scale dependence is also small. However, with decreasing £ it becomes moderate or even large, in 
particular for the phase. 
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1 Introduction 



Deeply virtual Compton scattering (DVCS) is considered as the theoretically cleanest process to 
investigate generalized parton distributions (GPDs) PQ El E]. These distributions are a hybrid 
of parton densities, form-factors, and distribution amplitudes and might be represented in terms 
of light-cone wave functions [U El E]- They are rather intricate functions depending on the 
longitudinal momentum fractions in the s- and t-channels, the momentum transfer squared, and 
the resolution scale. On the other hand they are phenomenologically very attractive, since they 
allow to combine information from different experiments in an optimal manner and since they 
encode non-perturbative information that cannot be extracted from either inclusive or elastic 
measurements alone. Their second moments provide e.g. the total angular momentum of partons 
in the nucleon [2] and the gravitational form factor of the nucleon. Moreover, a specific partial 
Fourier-transformation of GPDs is phenomenologically very interesting, as it provides functions 
which have a probabilistic interpretation. In the impact parameter space they can be viewed as 
parton densities in dependence of the longitudinal momentum fraction and the transverse distance 
from the proton center |H1 E] • The knowledge of transverse parton distribution does not only 
add substantially to our understanding of hadron structure but is also relevant for the prediction 
of cross sections in dependence of the impact parameter. For proton-proton scattering this has 
been especially emphasized with respect to LHC physics in Ref. |lUj . 

On one hand GPDs are thus a new window to study non-perturbative QCD and it has been 
already impressively demonstrated that they are experimentally accessible via the DVCS process 
[TT] IT!H ITU IT51 ITBj . On the other hand for several theoretical and experimental reasons 
the extraction of GPDs from measurements remains quite challenging because typically one is 
sensitive only to convolutions containing GPDs and one has to disentangle different contributions. 
The analysis simplifies substantially at large energies, where both photon EE EH] and vector- 
meson leptopro duct ion, e.g., Refs. [13 [TBI EH], have been measured by HI and ZEUS. In HERA 
kinematics the photon-proton interaction starts to be flavor blind and so one mainly accesses flavor 
singlet GPDs. Moreover, spin flip effects are suppressed, too, and only one set of GPDs is relevant, 
namely, the proton helicity conserved and parton helicity averaged GPDs H(x,£,t, Q 2 ). In this 
paper we concentrate on these singlet NNLO corrections, as the non-singlet case has already been 
studied in [20] ■ As usual, such analysis is only possible after adopting some parametrization for 
the dependence of GPDs on the s- and t-channel momentum fraction. We believe that realistic 
models can be most easily constructed by means of the partial wave decomposition of GPDs and 
amplitudes [2U E21 EH] , where the dominant contributions arise from the leading Regge trajectories 
[QUEUES- Support for this conjecture arises also from lattice calculations I27| I2H1 12§[ E?U|. 

In this letter we study radiative corrections to DVCS at and beyond next-to-leading (NLO) 
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order. This investigation is partially motivated by the fact that within a certain class of GPD 
models perturbative corrections at this order were reported to be rather large pH] 132] and so one 
should worry about the justification of the perturbative QCD approach. Recently, the radiative 
corrections in the flavor non-singlet case have been studied and it has been concluded that the rel- 
ative radiative corrections are moderate at NLO and become smaller at next-to-next-to-leading 
order (NNLO). Hence, these findings support the perturbative formalism in this sector. Consid- 
ering the singlet case at hand, we recall that the leading order (LO) contribution is given by the 
quark handbag diagram and that at next-to-leading order the gluon distribution appears as a 
new entry. It is known from deeply inelastic scattering (DIS) that the gluonic contribution in the 
small x-region is much larger than that of the sea quarks. Hence, the size of the NLO corrections 
depends in particular on the gluonic GPD and the appearance of large NLO corrections does not 
necessarily mean that perturbation theory fails. 

To clarify the situation, we employ here conformal symmetry to obtain the next-to-next-to- 
leading order corrections of the DVCS amplitude. In Sect.|2]we present, after a short introduction 
to the conformal approach, the analytic result for the DVCS amplitude in NNLO. In Sect. El relying 
on the pomeron pole as the dominant contribution at small momentum fraction, we numerically 
evaluate radiative corrections up to NNLO for the kinematics of HERA collider experiments. This 
analysis includes a comparison of the standard predictions in NLO with those of the conformal 
approach. We finally give our conclusions in Sect. 



2 The DVCS amplitude in NLO and NNLO 

The DVCS amplitude is defined in terms of the hadronic tensor 

T llv (q,P h P 2 ) = ^J d 4 xe i ^(P 2 ,S 2 \Tj,(x/2)j v (-x/2)\P 1 ,S 1 ), (1) 

where q = (q\ + q 2 )/2 (/j, and q 2 refers to the outgoing real photon). To leading power in Q 2 = —q\ 
(leading twist) and LO in the QCD coupling constant the hadronic tensor (PJ) is evaluated from 
the hand-bag diagram. In terms of the kinematical variables P = P\ + P 2 and A = P 2 — Pi, the 
result can be written as 

^ P ' A) = ~ -p—q-p—qj—q- ""'JpTtf + ' ' ' " ^ 

where the tilde-symbol denotes contraction = P MP X pcr V au with projectors 

V af3 = g af3 _ , ( 3 ) 

qi ■ q2 
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to ensure current conservation The ellipsis indicate terms that are finally power suppressed 
in the DVCS amplitude or are determined by the gluon transversity GPD, which is suppressed by 
a s /n and is not considered here. Note that to leading twist accuracy the parenthesis in (j2J can 
be replaced by g^. In the parity even sector the vector 

v° = u(p 2 , s 2 ) [nr + z^r) u{Pi, Sx) + ■ ■ ■ , (4) 

is decomposed into the target helicity conserving Compton form factor (CFF) H. and the helicity 
flip one S. Analogously, the axial- vector 

A° = U(P 2 , S 2 ) (W75 + Z^f) U(Pi, &) + •.., (5) 

is parametrized in terms of H and S, where again higher twist contributions are neglected. The 
normalization of the spinors is U(p, S) , y cr U(p, S) = 2p a . We also introduce the scaling variables 

P ■ q P ■ q 

where Q 2 = —q 2 . In DVCS kinematics and twist-two accuracy we have £ = r], while Q 2 = 2Q 2 . 

Before we proceed, let us decompose the CFFs, denoted by the set T = {TC, S, TC, £}, in flavor 
non-singlet (NS) and singlet (S) ones: 

F = Cfta™? + <f a *F , S ^ = ^ + G ^, (7) 

where the singlet piece contains the quark flavor singlet S JF and gluon G JF CFFs. The charge 
factors Q 2 with i = {NS, S} are given as linear combination of squared quark charges, e.g., the 
singlet one is given by the average of the squared charges for rif active quarks: 

Ql = ^Y,Ql (8) 

i=u,d,--- 

In the momentum fraction representation the Compton form factors are represented as convolution 
of the coefficient function with the corresponding GPD. In the singlet sector in which quark (^O) 
and gluon ( G t?) operators mix under renormalization, we might introduce the vector notation: 

f 1 dr 

A 2 , Q 2 ) = J i J C{x/t Q 2 //J?, a.(/i) \OF(x, V = A 2 , /i 2 ) . (9) 
Here the column vector 

G F 



F = {H,E,H,E} (10) 
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contains the GPDs, and the row one, defined as C = (^C, (l/£) G t Y ), consists of the hard scattering 
part that to LO accuracy reads 



We remark that the £ dependence in nC? and °C enters only via the ratio x/£. Note also that the 
w-channel contribution in the quark entry (jlljl has been reabsorbed into the symmetrized quark 
singlet distribution 



Here the second term in the square brackets with — (+)-sign for H, E (H, E)-type GPDs is 
for x > 7] related to the s-channel exchange of an anti-quark. The gluon GPDs have definite 
symmetry property under the exchange of x — ► — x: G H and E are even, while G H and E are 



The convolution formula (JJJJ) has already at LO the disadvantage that it contains a singularity 
at the cross-over point between the central region (—77 < x < 77) and the outer region (77 < x < 1), 
i.e., for x — £ = 77. Its treatment is defined by the ie prescription, coming from the Feynman 
propagator. The GPD is considered smooth at this point, but will generally not be holomorphic 
|34j . The fact that both regions are dual to each other, up to a so-called D-term contribution 
|35j . makes the numerical treatment even more complicated. This motivated our development 
of a more suitable formalism in The factorization scale \x in the GPDs is ambiguous and 

at LO order this induces the main uncertainty. Beyond LO this factorization scale dependence 
will be cancelled in the considered order of perturbation theory. The NLO corrections to the 
coefficient functions 36J and to the evolution kernels |3*7| EH] were predicted from conformal 
constraints, where the rotation to the standard MS scheme has been taken into account. Note 
that the conformal symmetry in MS scheme is broken and that the predicted results coincide with 
the diagrammatic evaluation |4"U1 141[ 14*2*] . To this order and in this scheme a numerical code has 
been made accessible that includes evolution, see, e.g., [32|- As already mentioned above, it was 
found that the perturbative corrections to NLO can be quite large. 

At present it seems hardly possible to study perturbative corrections beyond NLO accuracy 
in the standard scheme, since the diagrammatical evaluation would require enormous effort. For- 
tunately, we can employ conformal symmetry to relate the perturbative corrections at NNLO to 
those for DIS [HU El El], where the NNLO corrections in the vector case has been completed 
by the substantial effort of Vogt, Moch and Vermaseren [3B]- From these calculations we get the 
normalization of the Wilson coefficients and anomalous dimensions. The conformal predictions 
arise from the application of the conformal operator product expansion and are valid as long as 




(11) 




(12) 



odd. 
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the twist-two operators behave covariantly under conformal transformation |%7| EBJ EH]- This 
is certainly true at tree level and it also can be ensured for vanishing /3-function in any order 
of perturbation theory within a special renormalization scheme [48 . To make contact with the 
conformal OPE, we expand the hard-scattering amplitude in terms of Gegenbauer polynomials 
with indices 3/2 and 5/2 for quarks and gluons, respectively, and introduce the conformal GPD 
moments, which formally leads to 

oo 

A 2 , Q 2 ) = 2 r^C^QV/A «.(/*)) Fjit, ^, A (13) 

3=0 

The expansion coefficients Cj can be calculated by the projection: 

CAQ /^^)) = r(3/2TO + 4) 

x\J\c{x,&i l ?,«.m=i){ ( ^' +3)[1 - x2] ^ /2 3[1 _^5/ 2 J(*)- (14) 

Note that we have here rescaled the integration variable with respect to £ and that the integral 
runs only over the rescaled central region. The conformal moments of the singlet GPDs are defined 

as 

Fir, A 2 u 2 ) = r (3/2)rq + l) l [\ ^ ( rjCf \(x\ 
*M,A tf A) 2^ + 3/2)2 1,^ ^ o mcS&A") 

Here j is an odd (even) non-negative integer for the (axial-)vector case. 

In the forward kinematics (A —>■ 0), our conventions are such that the helicity conserved GPDs 
coincide with the flavor singlet quark distribution and with x times the gluon distribution. Hence, 
for the moments we have agreement with the common Mellin-moments of parton densities, e.g., 
for the helicity averaged GPD: 



lim H(x,ri) = ( | (x) , q- = \im Hj(ri) = [ dx x J ( 
n-o \xG J 3 v^o J y 




(16) 

Unfortunately, the series (fT3*j) does not converge for DVCS kinematics, in particular not in the 
outer region, and one has to resum the OPE [^TJ |2U| or, equivalently, one can use a dispersion 
relation |5U1 15T] . The result for s 7i in terms of a Mellin-Barnes integral reads 

/7r ' Xl C^/v^a^H^A 2 ^ 2 ). (17) 



A 2 , Q 2 ) = I J° dj r j -' [i + tan (?f\ 
In the following we write the perturbative expansion as 



C^Q 2 ///, Q 2 //i 2 , a s (ii r )) 
2^' +1 r(j + 5/2) 



r(3/2)r( J + 3) 



Cf + (Q 2 //, 2 ) + ^-Cf\Q 2 /n 2 , Q 2 //i 2 ) + 0{al) 



2tt ] y ' (2tt) 



where corresponding to our conventions, the LO Wilson coefficients are normalized as 



Cf = (l,0), (19) 

and here we choose to distinguish the renormalization (/i r ) and factorization (/x) scales. 

Let us first give here the DVCS NLO corrections in the MS scheme. We restrict ourselves 
to the analysis of the kinematically dominant contribution, i.e., s 7i, and so we provide here 
only the results for the vector case. The conformal moments (fPfj) can be obtained from Refs. 
|3TH I4D1 I4*H W2\ . Using the representation of Ref . jHT] , the integrals which are needed are evaluated 
in a straightforward manner for integer conformal spin, see Appendix C of Ref. (52! for the quark 
entries. The analytic continuation to complex j leads to 



T cf\Q/^ 2 )=C, 



*si(i +j )-l + i: 4S £ +1 U 



2 2(j + 1)(j + 2) (j + l)>0" + 2) a 



M 



+ ^1*^,(20) 

G r d) ro// ^_ o„ T (4 + gj + j 2 ) + SM + 2)] + 2 + 3j + j 2 Sq yf /i 2 
C,- - -2n/Z> (l + j)(2 + j)(3 + j) 2 Q2 ' ( 21 ) 

where CV = 4/3 and Tp = 1/2. The entries of the anomalous dimension matrix read at LO: 

^ = - c '( 3+ yTWT^- 4Sl0 + 1) ) (22) 

sg (o) A rp 4 + 3j + j 2 
in — —WfTp- — — -, 23 
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g?v (i> » or 4 + 3 J+ / /oil 



where j3o = 2rif /3 — 11(7^/3, Ca = 3. In the MS scheme also the complete anomalous dimension 
matrix is known to two- loop accuracy |q"3*1 I54j . However, the conformal moments will mix with 
each other and the solution of the evolution equation has so far not been given in terms of a 
Mellin-Barnes integral. 

The advantage of the conformal symmetry is that it predicts the Wilson coefficients. However, 

the symmetry is only valid in a special conformal scheme (CS 1 ). In such a scheme the structure 

lr The treatment of the terms proportional to f3, that break conformal symmetry, is of course ambiguous. We 
employ here the so-called CS scheme in which the running of the coupling is implemented in the form of the 
conformal operator product expansion (COPE) that is valid for a hypothetical fixed point. In particular, conformal 
moments are multiplicatively renormalizable to NLO. For details see Refs. |52l I2()| . 
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of the Wilson coefficients up to NNLO is 



cf\Q 2 /fj, 2 ) 



2 



sf\Q 2 /n 2 ) 



^3 'J 



J 1 3 ^ 1 3 



4 2) (Q 2 /^) (0) / (0) 



(26) 



(27) 



/-<(l)//-)2/ 2m Q 2 . 1 (0) (0), 2 Q 2 
H'r ^ r 



where the so-called shift coefficients s^\Q 2 /fi 2 ) can be expressed in terms of harmonic sums 
s p( n ) = ELi 1 /^ as 



sf\Q 2 //i 2 ) = S 1 (j + 3/2)-S 1 (j + 2) + 21n(2)-ln^, 
sf{Q 2 /n 2 ) = (sf{Q 2 /^ 2 )) 2 -S 2 (j + 3/2) + S 2 (j + 2), 
and c « = ( s c« G c«) are the Wilson coefficients known from DIS. We have at LO c\ 



Q 2 



at NLO 



Sl(\ + j) + + 2) - - + / 2 5 l(j) r - S 2 (j + f) 

iv j j 2 iu ; 2 2(j + 1)0 + 2) u ; 



g(D o T (4 + 3j+i 2 )5 1 (i) + 2 + 3j + i s 
c ■ — —zrifi f — 



(l+j)(2 + j)(3 + j) 



(28) 
(29) 

(1,0), 

(30) 
(31) 



and at NNLO they are given by the Mellin moments of the DIS partonic structure functions 
|43[ I44j . To simplify their evaluation, we take for a fit, given in ^2], rather than the exact 
expression. 

The evolution of the singlet (integer) conformal moments in this CS scheme is governed by 



a fOfL(°) , « 2 (^)-,(D , ^> ( 2 ) + (a 4 " 
2tt 1j (2tt) 2T? (27r) 3T? 



i^,A 2 ,/i 2 ^ 



i-2 



(32) 



where the mixing matrix A? is not completely known. In the vector case the anomalous dimen- 



sions are known to NNLO [56]. In absence of the mixing term, the solution of the renormalization 
group equation Fj(£, A 2 , /i 2 ) = £j(ft, /JLo)Fj(£, A 2 , /1q) can be given using the path-ordered expo- 
nential evolution operator 



r /-m da' ) 
£j(j*,Ho) =Vexpi- / — 7j .(a s (//')) ^ . 

I JfMO h 1 J 



(33) 
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In the numerical analysis we will only resum the leading logarithms and expand the non-leading 
ones 



.,b=± 



S ab a Pj + 



a s{v) 06.(1) 



2tt 



b Af'{n, im>) +-^ ab Af'(fi,to) + 0(a 



00 



Here the projectors on the {+, — } modes are 



±1 



'A, A, V o 



3 1 ' 



(34) 



(35) 



where the eigenvalues of the LO anomalous dimension matrix are 



±Ai = 1 (^(0) + ^ 1 ^(0) _ GG(0) 



v 2 V ,J 1 ' J y 1 2 



A straightforward calculation leads to the matrix valued coefficients 



4EG 7 (°)GS y (0) 



(36) 



2/3 c 



(37) 



c=± 
X C P,- 



ah 



Wo 



a P, 



2ft 7 ' 7i 



2/? 7 ^ 7 ^ , 



XWo|2) a -P, 



/3f - /?2/?0 (0) /?1 (1) (2) 



-7 



4/3 2 " 

where the /x dependence is accumulated in the following functions: 



2A) J 



7 + 7 



b Pj, (38) 



«6 



n(3 + % - % 



1 - 



a s (n) 



il3 + a \ j - b \ j 
00 



(39) 



The expansion coefficients of the /3 function are defined as 

- — — j Po + 7-j no Pi + /„ x 3 P2 + <A«J, 

(7 47T (47TJ 2 (471") d 

A) = /3 1 = i / -102, /3 2 : 



(40) 



325 9 5033 

H 

54 f 18 



2857 



The expansion of the evolution operator (|34j) will then be consistently combined with the Wilson- 
coefficients (jTSjl . see for instance Ref. jS21- 
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3 Numerical evaluation of radiative corrections 



Here we study the perturbative corrections in the small-£ region. We adopt a simple ansatz for 
the conformal moments that is inspired by the dominance of the pomeron and by the assumption 
that 0(£ 2 ) terms in the conformal polynomials are insignificant: 

H (£ A 2 Q 2 ) = ( Ar s ea sea F(A 2 )B(l+j-a sca (A 2 ),8)/B(2-« sca (0),8) \ 
' \ iV G G F(A 2 )B(l+j-a G (A 2 ),6)/B(2-a G (0),6) J 

where B(x,n) = T(x)T(n)/T(x + n), and the ellipsis denotes the neglected 0{£ 2 ) terms, as well 
as valence components whose contributions are also small for small £. Here the normalization 
is sea F(A 2 = 0) = G F(A 2 = 0) = 1 and for aj(A 2 ) we use the "effective" pomeron trajectory 
a (A 2 ) = a(0) + a' A 2 . We remind that in deeply inelastic scattering the structure function 
F2 ~ (l/xBj) A( '^ 2 ' ) grows with increasing Q 2 . Here the exponent is related to the intercept of the 
Regge trajectory A = a(0) — 1. The values for a(0) will be specified below. Although also the 
slope a 1 is scale dependent [211, we choose here the standard value of the soft pomeron a' = 0.25. 
In the forward case the moments (J41j) arise from the parton densities 

A ' ' ■ x- a ^ (1 - x) 1 + uy(x) + d v (x) , G = — Na . n . -r x~ a ^ (1 - a;) 5 , 



B(2-a sea (0),8) v ' vw vw ' B(2-a G (0),6) 

(42) 

for which we have adopted a generic realistic parametrization. Note that the valence component 
u v + d v is not taken into account in Eq. (|41j) . since it is a non-leading contribution for small x. 
The normalization factors are related by the momentum sum rule 

I dxx[E(x) +G(x)} = 1 N G + N sca + f dx x[u v + d v ](x) = 1 . (43) 

Jo Jo 

In the asymptotic limit Q — ► 00, the evolution equation tells us that iV G = 4Cf/(4:Cf + n/) i.e. 
that more than 50% of the longitudinal proton momentum is carried by gluons. At experimentally 
accessible large scales the gluons already carry about 40 % of the momentum. For the momentum 
of the valence quarks we choose the generic value 1/3 and so A^ sea = 2/3 — N G . 

Usually, the GPDs are taken from some non-perturbative (model) calculation or ansatz and 
plugged into the CFFs at a given input scale. Since the scale and the perturbative scheme are 
usually not specified, the matching of perturbative and non-perturbative frameworks has its own 
uncertainties. Let us first study the changes of the CFF (fTTjl in a given scheme and input scale that 
appear when one includes the next order. The changes to the modulus and phase are appropriately 
measured by the X-factors: 

,N p LO\ 



In 

K = - 
In 



S^N^-iLO 



arg ( S H^ 



fCf = ^ ^ . (44) 




Figure 1: The relative radiative corrections, defined in Eq. (I44j) . are plotted versus £ for the 
logarithm of the modulus [(a) and (b)] and phase [(c) and (d)] of S TC, see Eqs. (jT7jl and (|41|). for 
A 2 = [(a) and (c)] and A 2 = -0.5 GeV 2 [(b) and (d)]: NNLO (solid) as well as in NLO for 
the CS (dashed) and MS (dotted) scheme. Thick (thin) lines refer to the "hard" ("soft") gluon 
parameterization and we always set ji — Q and a s (Q 2 = 2.5 GeV 2 ) /n = 0.1. 

Here S 7^ N LO denotes the N p LO approximation, e.g., P = for LO. One should bear in mind 
that fT-factors actually measure the necessary reparameterization of the GPD to fit the given 
experimental data. We take the ansatz (|4"T|) with sca F(A 2 ) = G .F(A 2 ) and so the factorized A 2 
dependence essentially drops out in the .ff -factors. The ratio of gluon GPD to quark one is 
controlled by the factor NQ/N seau and, more importantly, by the differences of intercepts O£q(0) — 
a sea (0). To study the influence of this ratio, we distinguish two cases: 

H) "hard" gluon: N G = 0.4, a G (0) = a sea (0) + 0.1 , (45) 

S) "soft" gluon: N G = 0.3, a G (0) = a sca (0) . (46) 

We will use these parameters and a sea (0) = 1.1 at the input scale Q 2 = 2.5 GeV 2 . Moreover, we set 
[i — Q and independently of the considered approximation we choose a s (/z 2 = 2.5 GeV 2 ) = 0.17T 
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and set the number of active flavors to three. 

In Fig. Q we depict for the typical kinematics of HERA collider experiments, i.e., 10~ 5 < £ < 
5 ■ 1CT 2 , the resulting K factors for the logarithm of the modulus [(a) for A 2 = and (b) for 
A 2 = -0.5 GeV 2 ] and phase [(c) for A 2 = and (d) for A 2 = -0.5 GeV 2 ]. Here the thick and 
thin lines correspond to the "hard" and "soft" gluon parameterizations, respectively. We observe 
an almost flat £ dependence of the K\ factors in panels (a) and (b). This is not surprising, since 
the essential contribution arises from the pomeron pole and the CFF behaves as: 

S n~(-J i + tan(|(a(A 2 ) - 1)) In \ S H\ ^ «(A 2 ) ln(l/f) + const. . (47) 

For small £ this leads to the flatness we observe. The size of perturbative NLO corrections, see 
dashed (CS scheme) and dotted (MS scheme) lines, essentially depends on the ratio of gluon 
to quark GPDs. Since the gluons are a new entry, formally counted as NLO contribution, this 
finding is obvious and goes along with the observation that the perturbative corrections strongly 
vary within the used parameterization of parton densities in the Radyushkin GPD ansatz [35] . In 
the "hard" gluon parameterization the logarithm of the modulus reduces about 7-11% [5-8%] in 
the MS [CS] scheme, corresponding to the reduction of the modulus itself in the range of 40-70% 
[30-55%], where the drastic upper values correspond to £ = 10~ 5 . 

The relative radiative corrections to the phase grow in the small £ region with decreasing £ 
and can be of the order of up to 24% [13%] in the MS [CS] scheme. These effects are related 
to the signs for NLO Wilson coefficients, see Eqs. (|2D jl . (j2T j) . (jSIT j) . and (ETTjl. In the "soft" gluon 
parameterization the NLO corrections are quite moderate for the modulus [(a) and (b)] and 
negligible for the phase [(c) and (d)]. From all four panels it can be realized that compared to MS 
scheme in the CS one the NLO corrections are typically reduced by 30-50%. This reduction has 
been also observed in the flavor non-singlet case [20] ■ The NNLO corrections (solid), compared to 
the NLO (dashed) ones, are drastically reduced. For the "soft" gluon parameterization they are 
practically negligible while for the "hard" gluon input they are reduced to the 1-2% level, except 
for the phase with A 2 = — 0.5GeV 2 , where 5% are reached at £ = 10~ 5 . 

Let us finally address the modification of the scale dependence due to the higher order cor- 
rections. We only consider here the CS scheme and, analogously as in Eq. we quantify the 
relative changes due to the evolution by the ratios 

/NP-iLO' 



<iln 



S H NPL ° /din 



darg(V PLO ) /daig^H 



A d\nQ 2 I c/lnQ 2 ' ^ arg dlnQ 2 / din Q 2 ' ^ 

For the (exact) evolution of a s (Q) we take the same scale setting and initial condition as above. 
However, the conformal moments (|41j) are evolved in the CS scheme, starting at the input scale 



11 




Figure 2: The relative change of scale dependence, cf. Eq. (jlHj) . in the CS scheme at NLO (dashed, 
dotted) and NNLO (solid, dash-dotted) versus £ is depicted for the logarithm of the modulus (a) 
and phase (b) of the CFF (JUJ) with A 2 = (dashed, solid) and A 2 = -0.5 GeV 2 (dotted, dash- 
dotted) and Q 2 = 4 GeV 2 . We set \i = Q, a s (/i 2 = 2.5GeV 2 )/7r = 0.1 and took the input flU} at 
the scale Q 2 , = 1 GeV 2 . Thick and thin lines correspond again to "hard" and "soft" gluonic input. 

Qq = 1 GeV 2 , to Q 2 = 4 GeV 2 . The non-leading logs in the solution of the evolution equation 
are expanded with respect to a s and are consistently combined with the Wilson-coefficients 
()18|) in the considered order. The unknown NNLO mixing term A^ fc in Eq. ()32|) is neglected. This 
mixing can be suppressed at the input scale by an appropriate initial condition and so we expect 
only a minor numerical effect; see also Ref. [HE]- The dashed and dotted lines in Fig. |2] show that in 
NLO the scale dependence changes can be rather large even of about 100% or more. In general the 
relative radiative corrections to NNLO are getting smaller. For instance, the NNLO corrections in 
panel (a) are almost negligible for the "soft" gluonic input with A 2 = (thin solid), they increase, 
however, for A 2 = —0.5 GeV 2 and are becoming large for the "hard" gluonic input, e.g., about 
—35% at £ = 10~ 5 (thick dash-dotted). Note that these large corrections at very small £ are 
essentially caused by those of the anomalous dimensions in the vicinity of j = 0, corresponding to 
the large corrections of the gluon splitting kernels at small x, reported in The same sources 
also cause the huge NNLO corrections to the phase in panel (b). We remark that the modulus 
of s 7i is dominated by its imaginary part for which radiative corrections are milder than for the 
real part. The real part and so also the phase at very small £ are rather strongly affected by the 
NNLO corrections to anomalous dimensions. On the other hand for 5 ■ 10~ 4 < £ and 5 ■ 10~ 3 < £ 
the radiative NNLO corrections to the logarithm of the modulus (a) and phase (b), respectively, 
are rather mild (solid and dash-dotted lines). Restricted to these kinematics our findings support 
the convergence of the perturbative series. 
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4 Summary 



In this letter we have studied NLO and NNLO corrections to deeply virtual Compton scattering 
in the small £ region. We confirmed that large radiative corrections at NLO can appear, reported 
before, and clarified their source which is entirely tied to the gluonic sector. In particular, if 
the gluon distribution starts to have a steeper increase at small £ than the quark ones, the NLO 
corrections will be dominated by the negative NLO gluon contribution and so the modulus of 7i 
will drastically reduce. On the other hand, if the gluon contribution is relatively small, already 
the NLO corrections are moderate. In any case the NNLO corrections are becoming moderate 
or even small at a given input scale, even at a few GeV 2 . This fact supports the perturbative 
framework of DVCS. 

The situation with respect to the scale dependence is not so conclusive. Going from LO to NLO 
we observe in general a big enhancement that arises from the large corrections to the anomalous 
dimensions, cf. |46| . To NNLO they will be reduced and the relative changes for the logarithm 
of the modulus are getting reasonable but grows to be large with decreasing £. Note that in this 
region the NNLO gluonic evolution effects are comparable in size with the NLO ones Also 
the NLO radiative corrections to the scale dependence of the phase of s 7i are rather large, in 
particular for A 2 = 0, at the scale of 4 GeV 2 . To NNLO accuracy the convergency improves for 
5 • 1CU 3 < £. Unfortunately, at smaller values of £ the convergency is lost. These large corrections 
due to evolution at small £ are certainly related to those found in DIS j4Tj . 

If one is interested to access GPDs from the DVCS cross section measurement at small £, only 
the modulus of s 7i is essential. In that case perturbation theory seems to work in the sense that 
NNLO corrections of the Wilson-coefficients are negligible. They are, however, important for the 
scale violating effects for £ < 5 • 10 -4 (at relatively low Q 2 ~ 4 GeV 2 ). We also conclude that 
the photon and vector-meson leptoproduction data taken by the HI and ZEUS collaborations 
should be perturbatively analyzed at NLO |57|. To our best knowledge a common perturbative 
analysis has not been done so far. To achieve this in a simple and numerical stable manner, the 
Mellin-Barnes integral representation seems to be preferred. 

This project has been supported by the U.S. National Science Foundation under grant no. 
PHY-0456520, German Research Foundation (DFG), and Croatian Ministry of Science, Education 
and Sport under the contract no. 0119261, as well as The National Foundation for Science, Higher 
Education and Technological Development of the Republic of Croatia under the contract 01.03./02. 
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